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Weibull 分 布 兴趣 参数 的 广义 置信 区 间 * 


(北方 工业 大 学 理学 院 统计 系 ， 北 京 100144) 

摘 要 : Weibull 分 布 可 以 用 来 描述 疲劳 失效 、 真 空 管 失效 和 轴承 失效 等 ， 是 重要 的 寿命 分 布 。 本 文 研究 
了 Weibull 分 布 尺 度 参数 、 形 状 参数 、 分 位 数 及 可 靠 度 函数 的 广义 置信 区 间 问 题 。 利 用 广义 枢 轴 
量 给 出 四 个 兴趣 参数 的 广义 置信 区 间 。 证 明了 由 广义 枢 轴 量 确定 的 四 个 兴趣 参数 的 广义 置信 区 间 
具有 频率 意义 下 的 实际 置信 水 平 。 
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Weibul 分 布 是 应 用 广泛 的 寿命 分 布 之 一 。 对 于 Weibul 分 布 的 兴趣 参数 的 置信 区 间 ， 前 人 
做 过 大 量 的 研究 也 提出 了 许多 方法 。 诸 如 基于 参数 的 最 好 线性 无 偏 估 计 、 最 好 线性 同 变 估计 或 
最 大 似 然 估计 确定 未 知 参 数 的 精确 置信 区 间 ， 用 条 件 分 布 的 方法 给 出 参数 的 置信 区 间 钙 ， 这 些 
方法 在 计算 上 非常 复杂 ， 不 便于 处 理 。Rekkasa 和 Wongb 运用 似 然 分 析 理 论 给 出 参数 的 近似 置 
信 区 间 四 ， 这 种 方法 虽然 有 较 高 的 精确 度 ， 但 是 计算 仍然 很 复杂 。 本 文 运用 Weerahandi 提出 
的 广义 枢 轴 量 方法 回 给 出 Weibull 分 布 尺度 参数 、 形 状 参数 、 分 位 数 和 可 人 靠 度 函数 的 精确 置信 
区 间 。 

当 我 们 求 参 数 置 信 区 间 时 ， 如 果 有 讨厌 参数 存在 ， 获 得 参数 的 精确 置信 区 间 将 存在 困难 ， 
对 此 Weerahandi 提 出 了 广义 枢 轴 量 的 概念 用 于 求 参数 的 置信 区 间 , 令 R=7r(X;z,é&) 是 X,z,é 
的 函数 ， 其 中 上 = (0,7)， 若 及 满足 以 下 两 条 性 质 : 

1) RR 的 分 布 与 & = (67 无关; 

2) 及 的 观测 值 robs = 7(zi x,é) 不 依赖 于 讨 大 参数; 

则 称 RR 为 广义 枢 轴 量 ， 如 果 RR 的 样本 空间 子 集 C; 满足 等 式 P(Re Cr) = ?7(0 < YY < 1)， 则 参 
数 空 间 的 子 集 9G.(7r) = {0 € ol ZiE) e Cr} 称 为 6 的 1007 名 的 广义 置信 区 间 。 

本 文 结构 如 下 第 二 节 给 出 兴趣 参数 的 广义 置信 区 间 ， 第 三 节 证 明 由 广义 枢 轴 量 确定 的 广 

义 置信 区 间 具 有 频率 性 质 。 


2 Weibull 分 布 兴 趣 参 数 的 广义 置信 区 间 
随机 变量 Y 服从 Weibull 分 布 ， 概 率 密度 为 
B/y\o-!i VAN 
RS 骂人 站 Y>0, 
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其 中 有 > 0 是 形状 参数 ， 和 > 0 是 尺度 参数 。 令 X= InY, 入 = en， B= 则 


Xem{f se}om{-es}=3/(S*) 
这 里 f(t) = ete-e， 一 00<t< +o0。 所 以 X 服 从 极 值 分 布 ， 是 位 置 -尺度 分 布 族 。 
设 Xi, 2,… ,Xn 是 来 自 此 总 体 的 独立 同 分 布 样本 。 令 


n 


-， 本 二 
六 = > MX，92= >》 (Xi 一 部)2. 


?一 1 ?一 1 


徐 兴 忠 、 李 国 英 在 枢 轴 分 布 族 中 的 Fiducial 推断 外 给 出 广义 枢 轴 模 型 况 芋 /poB， S$ ob, 


其 中 n nn 
Bi =i De E2 = 》 (ei— Ei)’, ei 7 一 上 2 ,nN 
Se 一 
是 来 自 标 准 极 值 分 布 的 独立 同 分 布 样本 。 则 Fiducial 模 型 为 = 五 -BE o = 高 。 据 此 给 
出 Weibull 分 布 尺度 参数 、 形 状 参数 、 分 位 数 和 可 靠 度 函 数 的 广义 枢 轴 量 。 
考虑 Weibull 分 布 尺 度 参数 入 的 100(1 - oa) 名 广义 置信 区 间 (0 < a < 1)。 构 造 广义 枢 轴 量 


Ri ="{(X,8);(2,s), (1,0)} =£— 寻 ， 


其 中 忆 ; = 2>k，Ezs = 3， 元 s 是 义 , 5 的 观测 值 。Ri 的 分 布 与 参数 jo 无关 ， 并 且 7{(5,s); 
(元 , s), (1, Ea = jy， 满 足 广 义 枢 轴 量 的 两 个 条 件 。 设 0 < a < 1, 令 抽 7, 910 分 别 表 示 Ri 的 分 
布 的 1009 免 分 位 数 和 100(1 一 全 )% 分 位 数 ,有 P{91L < Ri <jv}=1 一 a; 则 ( 抽 4,91v) 是 jy 的 
100(1 一 a)% 广 义 置 信 区 间 。 所 以 Weibull 分 布 尺 度 参数 的 100(1 一 a)% 广 义 置信 区 间 为 
( exp{$15}, exp{f1v}). (1) 
考虑 Weibull 分 布 形状 参数 8 的 100(1 一 a)%% 广 义 置信 区 间 。 构 造 广义 枢 轴 量 Ro = 7{5; s, B} 
= 人 笃 ， 其 中 机 = 2 = S86。 显然 Ra 满足 广义 枢 轴 量 的 两 个 条 件 。 令 9oL, 92v 分 别 表 
示 Rs 的 分 布 的 1009 久 分 位 数 和 100(1 一 )%% 分 位 数 ， 有 Pf{9oL < R2 < 92v} =1-a 
则 8 的 100(1 一 a)% 广 义 置信 区 间 为 
(G2L, G2U). (2) 
考虑 Weibull 分 布 分 位 数 的 100(1 一 a)% 广 义 置信 区 间 。Weibull 分 布 的 分 位 数 为 [一 ln(1 一 
Qa)]， 可 以 写 为 ex[ 一 In(1 -~ ojj"。 构 造 广义 枢 轴 量 


Rs =7{ (X, 5); (3,3), (100) } =e’ 副 [In(1 -a)] 记 ， 


Ra 满足 广义 枢 轴 量 的 两 个 条 件 。 令 83, ga 分 别 表 示 Rs 的 分 布 的 1008 色 分 位 数 和 100(1 一 
9)% 分 位 数 ， 有 P{5sL < Rs < 93v} = 工 - a， 则 分 位 数 和 [一 In(l 一 a)]# 的 100(1 一 a)% 广 义 
置信 区 间 为 
(G3L, 93U). (3) 
考虑 Weibull 分 布 可 靠 度 函数 的 100(1 一 a)%% 广 义 置信 区 间 。 在 给 定 yo 时 ，Weibull 分 布 可 
靠 度 函数 为 exp{ 一 ( 实 )8}， 令 zo = ljnyo， 则 可 靠 度 函 数 为 exp{-e s}， 不 妨 先 求 加 =& 的 广 
义 置 信 区 间 。 构 造 广义 枢 轴 量 


Ra = r{ (X, 5); (2, s), (1,0)} 3 2 
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可 以 证 明 满足 广义 枢 轴 量 的 两 个 条 件 。 令 8941, 94v 分 别 表示 Rs 的 分 布 的 1009% 分 位 数 
和 100(1 一 全)% 分 位 数 ， 有 P{94r < Ra < 9av} = 1 一 Q， 则 (94z, 940) 就 是 开 兴 的 100(1 一 
a)% 广 义 置信 区 间 。 所 以 可 靠 度 函数 exp{e”s} 的 100(1 一 a)% 广 义 置信 区 间 为 


( exp{—e30 }, exp{—edr} ). (4) 


注 1 徐 兴 忠 ， 李 国 英 向 运用 Fiducial 推断 的 方法 也 给 出 了 可 靠 度 函 数 的 精确 置信 限 ， 其 结 
果 与 本 文 据 广义 枢 轴 量 确定 的 可 靠 度 函数 的 广义 置信 区 间 是 完全 一 致 的 。 

注 2 本 文 给 定 的 Weibul 分 布 兴趣 参数 的 广义 置信 区 间 可 以 用 Monte Carlo 方 法 计算 得 
出 ， 因 此 与 前 人 的 结果 相 比 ， 运 算 过 程 更 加 简捷 方便 。 


3 ”频率 性 质 


本 节 研 究 广义 置信 区 间 的 频率 性 质 就 是 研究 兴趣 参数 落 在 广义 置信 区 间 的 覆盖 概率 。 

定理 1 Weibul 分布 尺 度 参 数 入 的 广义 置信 区 闻 (1) 具有 频率 意义 下 的 实际 置信 水 平 1 一 
a， 即 P{exp{gir} < 入 < exp{giu}} =1-a。 

证 明 


P{ exp{ffz} < 入 < exp{gir}j} = Prz <A<Dol=Pf<brl-Pfw<9zrh， 


P{p < 91v} = P{ER(o < 1-3}= P{ P(X -5 <4) <1-} 


这 里 FR, (.) 表示 Ri 的 分 布 函 数 。 设 羡 =j 十 oBY, 5S = o， 其 中 忆 +, 包 分 别 与 Bi, Bz 独立 
同 分布 。 令 Q* 表示 (B+, 到?) 的 分 布 函 数 ，Q@ 表示 (Bi, E2) 的 分 布 函数 。 则 上 式 为 


cfo( 时 = 熏 } < 中 -afe( 品 )> 外 -1 和 


这 里 下 (.) 表 示 盈 的 分 布 函数 。 所 以 P{exp{giz} < 和 < exp{fpiz}} = 1 一 Qa， 结论 得 证 。 

定理 2 Weibul 分 布 形状 参数 6 的 广义 置信 区 间 (2) 具有 频率 意义 下 的 实际 置信 水 平 1 一 
a， 即 P{bar <B< Ov}=1—oa。 

证 明 证 明 同 定理 1， 故 略 。 

定理 3 Weibull 分 布 分 位 数 和 [一 ln(1 - a)]# 的 广义 置信 区 间 (3) 具有 频率 意义 下 的 实际 置 
信 水 平 1 -a, 即 Pf{63r < 和-In(1 -oa)]3 < 93v}=1-a。 

证 明 


P{ $3 < MN- In(1 — o)]# < Gav } 


= P{A~In(1 — oa)]? < 930} — P{A-In(l ~ o)]# < $3 }. 
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与 定理 1 做 相同 的 标记 ， 则 
P{A—In(1 -oa 二 < sv } = P{F(A-In(1 —a)]$)<1- 2) 


P{P(ex 5 刍 [-In(1 -a)] 训 <A[- ni-o)]?)<1-?} 


| 


Q"{o{ < El 0 上 i 


E? — In[— ln(1 — Q)] Ql a 
GT 索 J 
这 里 F() 表示 全 于 愉 二 9 的 分 布 函数 ， 所 以 P{93L < 和 -In(1 -oa 了 <93v}=1-a， 结 
论 得 证 。 

定理 4 Weibull 分 布 可 靠 度 函数 exp{ 一 (各)3} 的 广义 置信 区 间 (4) 具有 频率 意义 下 的 实际 
置信 水 平 1 一 a， 即 


P{ exp{—e?+} < exp { 一 ( 半 | < exp{ -e307}} =1—a. 


证 明 证 明 同 定 理 1， 故 略 。 

由 上 述 四 个 定理 可 知 广义 枢 轴 量 确定 的 Weibull 分 布 兴趣 参数 的 广义 置信 区 间 的 覆盖 概率 
为 1 一 a。Rekkasa 和 WongbDl 由 似 然 比方 法 给 出 了 兴趣 参数 的 近似 置信 区 间 不 具备 这 样 的 性 
质 。 


参考 文献 : 


[1] Lawless J F. Statistical Models and Methods for Lifetime Data[lM]. New York: Wiley, 1982 

[2] Rekkas M, Wong A. Third-order inference for the Weibull distribution[Jj. Computational Statistics & Data 
Analysis, 2005, 49(2): 499-525 

[3] Weerahandi S. Generalized confidence intervals[J]. Journal of the American Statistical Associaation, 1993, 
88(423): 899-905 

[4j 徐 兴 忠 ， 李 国 英 . 枢 轴 分 布 族 中 的 Fiducial 推 断 [J]. 中 国 科学 A 辑 ， 数 学 ，2006, 36(3): 340-360 
Xu X 2, Li G Y. Fiducial inference in the pivotal family of distributions[J]. Science in China: Series A， 
Mathematics, 2006, 49(3): 410-432 


Generalized Confidence Intervals for Interest Parameters of the 
Weibull Distribution 
ZHAO Gui-mei, CUI Yu-jie 


(College of Science, North china University of Technology, Beijing 100144) 

Abstract: The Weibull distribution, a kind of important life distributions, can be used to describe 
fatigue failure, vacuum tube failure and bearing failure, etc. In this paper, the generalized confidence 
intervals for the scale parameter, the shape parameter, the percentile and the reliability function of the 
Weibull distribution are investigated. Generalized confidence intervals of four interest parameters are 
constructed by using the concept of generalized pivotal quantity. At last, we prove that they have the 
exact confidence levels in sense of frequency. 
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